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A  POSSIBLE  EXPERIMENTAL  DESIGN 
FOR  TESTED  THE  EFFECT  OF  ARTIFICIAL  CLOUD  SEEDING 


Because  of  the  inherent  variab'  ity  of  rainfall,  it  is  very  difficult 
to  determine  ■whether  or  not  the  artificial  seeding  of  clouds  has  produced  pre¬ 
cipitation  over  and  above  that  which  would  be  expected  naturally.  Current  ar¬ 
guments  pro  and  con  attest  to  this  difficulty.  Much  disagreement  can  be  avoid¬ 
ed  if  the  rules  of  experimental  procedure  and  interpretation  are  laid  down  in 
advance,  for  controversies  are  more  difficult  to  resolve  if  the  data  are  analyz¬ 
ed  after  the  fact.  For  this  reason,  we  have  endeavored  to  develop  an  experimen¬ 
tal  design  of  such  nature  that  a  given  outcome  can  be  appraised  objectively  on 
the  basis  of  the  probability  of  obtaining  such  a  result  by  chance,  if  seeding 
were  assumed  ineffective. 

For  purposes  of  this  report,  we  consider  it  unnecessary  to  present  de¬ 
tailed  masses  of  meteorological  information,  since  much  of  this  is  common  know¬ 
ledge  to  all  who  have  analyzed  weather  records.  Drawing  upon  well  established 
facts,  wc  propose  to  discuss  the  logical  and  statistical  issues  involved  in  the 
design  of  experiments  capable  of  determining,  with  reasonable  assurance,  whether 
or  not  seeding  causes  a  significant  increase  in  rainfall.  Besides  giving  expli¬ 
cit  directions  for  applying  the  technique  which  we  have  devised,  we  cite  evi¬ 
dence  that  its  underlying  assumptions  are  sound  and  show  that  it  is  sensitive 
enough  to  detect  even  a  relatively  small  increase  over  natural  precipitation* 

If  it  were  possible  to  predict  with  perfect  accuracy  how  much  preci¬ 
pitation  would  occur  in  the  absence  of  seeding,  it  would  be  a  simple  matter  to 
determine  the  effect  of  seeding.  Again,  even  without  an  accurate  forecast  of 
rainfall  aside  from  artificial  modifications,  it  might  still  be  possible  to  e- 
valuate  the  contribution  due  to  seeding,  if  the  ensuing  physical  processes  were 
procisely  known.  Unfortunately,  neither  condition  is  satisfied  at  present,  and 
in  default  of  both,  the  problem  of  testing  the  effect  of  seeding  is  statistical. 
The  natural  pattern  of  rainfall  is  extremely  variable  from  one  storm  to  another, 
and  even  monthly  averages  are  highly  divergent  from  year  to  year;  nevertheless, 
under  natural  conditions,  there  exists  a  certain  probability  distribution  of 
rainfall.  If  seeding  wore  effective,  the  observed  rainfall  would  show  a  greater 
tendency  toward  larger  values  than  under  natural  conditions.  The  technique  of 
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deterndning  whether  or  not  a  real  change  has  takpn  place  cones  under  the  general 
heading  of  the  testing  of  statistical  hypotheses.  The  province  of  experimental 
design  is  to  provide  settings  wherein  valid  tests  can  be  made. 

The  statistical  theory  of  experimental  design  boils  down  to  the  con¬ 
struction  of  a  test  variato  satisfying  two  conditions  s 

1,  The  test  variate  measures  the  relevant  effect  of  the  treatment  in  ques¬ 
tion. 

2.  Its  probability  distribution  is  known  under  the  null  hypothesis  that 
the  experimental  treatment  docs  not  affect  the  test  variate. 

In  ordinary  analyses  of  observations  drawn  from  a  specified  population,  the  task 
of  constructing  an  appropriate  test  variate  is  relatively  simple,  since  one  can 
usually  hypothesize  the  independence  of  individual  observations,  or  at  wor3t,  a 
form  of  dependence  which  presents  no  analytical  obstacles.  The  aim  of  the  tech¬ 
nician  is  then  to  choose  from  several  valid  tests  one  which  has  maximum  sensi¬ 
tivity.  Having  performed  the  experiment  and  computed  the  statistic  previously 
decided  upon,  one  passes  judgement  as  to  the  tenability  of  the  null  hypothesis. 
If  the  value  of  the  to3t  variate  i3  found  to  lie  very  far  outside  of  the  corn- 
iron  interval  of  variation,  the  null  hypothesis  iB  rejected,  and  one  concludes 
that  a  significant  effect  has  taken  place. 

Let  us  review  the  two  conditions  upon  which  a  valid  test  depends. 

First,  the  measurement  must  be  relevant.  This  requirement  rules  out  both  in¬ 
cidental  phenomena  and  significant  effects  apart  from  the  experimental  treatment. 
Any  action  whatever  must  have  some  kind  of  result,  but  unless  the  result  obtain¬ 
ed  is  of  the  kind  intended,  it  is  beside  the  point.  For  exainplo,  when  13  play¬ 
ing  cards  arc  dealt  from  a  standard  deck  without  jokers,  it  i3  certain  that  the 
outcome  will  be  a  legitimate  hand  of  bridge.  Now,  the.  probability  of  obtaining 
any  preassigned  hand  as  a  random  draw  from  tho  deck  is  about  10“ 22,  and  if  the 
hand  were  specified  in  advance,  tho  event  of  drawing  it  would  be  overwhelming¬ 
ly  significant  of  skill.  Unforsecn  as  a  particular  entity,  the  same  hand  sig¬ 
nifies  nothing  more  nor  less  than  the  fact  that  some  permutation  is  bound  to 
occur.  A  continuous  variate  presents  this  situation  in  its  logical  extreme. 
Here,  provided  the  distribution  function  is  continuous,  tho  probability  of  ob¬ 
taining  ary  pro assigned  value  exactly  is  aero*  therefore,  ono  considers  the 
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probability  of  obtaining  a  value  v/ithin  a  stated  neighborhood  of  a  given  point, 
Without  roforoncc  to  actual  data,  it  must  bo  decided  how  the  test  variate  will 
be  interpreted  with  regard  to  the  intended  operation  of  the  treatment.  This  is 
a  logical  issue  which  can  be  settled  a  priori.  For  instance,  if  the  purpose  of 
the  treatment  is  to  increase  a  certain  variable,  then  it  is  appropriate  to  in¬ 
terpret  the  experimental  result  in  relation  to  the  class  of  all  values  equal  to 
cr  greater  than  the  one  obtained.  Once  an  admissible  standard  of  interpretation 
has  been  defined,  and  a  given  experimental  outcome  proves  to  be  highly  improb¬ 
able  on  the  basis  of  the  null  hypothesis,  the  rejection  of  the  null  hypothesis 
is  then  reasonable.  On  the  other  hand,  a  significant  influence  might  be  exert¬ 
ed  by  some  dominant  factor  other  than  the  treatment  deliberately  applied. 

■Where  it  is  felt  that  a  situation  of  this  sort  might  exist,  great  care  must  be 
exercised  in  the  construction  of  the  test  variate,  in  order  that  the  treasured 
effect  be  properly  attributable  to  the  experimental  treatment  alone.  This  is 
one  of  the  most  vulnerable  points  in  the  design  of  experiments.  Techniques  of 
controlling  or  at  least  accounting  for  extraneous  sources  of  variability  are  in¬ 
troduced  to  prevent  inflation  of  the  supposed  effect  of  the  experimental  treat¬ 
ment.  Such  techniques,  when  correctly  employed,  have  tho  further  advantage  of 
making  the  test  more  sensitive.  In  addition  to  direct  methods  of  control,  a 
safeguard  against  unsuspected  or  unmanageable  factors  is  supplied  by  randomiz¬ 
ation.  Randomization  does  not  eliminate  disturbing  influences  but  corrects  for 
them  in  a  probability  sense  try  affording  an  opportunity  for  the  test  variate  to 
be  decreased  as  often  as  increased  by  their  presence.  A  more  cogent  reason  for 
randomization,  however,  is  that  tho  whole  structure  of  the  statistical  testing 
of  hypotheses  is  built  upon  the  concept  of  random  sailing. 

The  second  condition  for  validity  is  that  the  probability  distribution 
of  the  test  variate  be  known  under  the  null  hypothesis .  To  meet  this  require¬ 
ment,  it  is  ordinarily  necessary  to  proceed  in  steps.  The  measurements  that  can 
be  made  directly  are  usually  found  to  involve  several  sources  of  variation  at 
once.  The  first  step,  therefore,  is  to  design  the  experiment  in  such  a  way  that 
the  net  effect  of  the  experimental  treatment  can  be  derived  from  the  direct  mea¬ 
surements  .  As  a  rule,  tho  probability  distribution  of  the  synthetic  variable 
thus  constructed  will  be  known  theoretically  in  its  general  character,  but  there 
will  be  some  unknown  parameters  which  will  have  to  bo  specified  before  the  dis¬ 
tribution  function  can  be  evaluated  quantitatively.  Almost  without  exception, 
however,  it  is  impossible  to  determine  the  exact  value  of  a  paramo  tor  from  ob- 


servational  data.  Therefore,  the  statistician  proceeds  to  the  next  stage  of 
his  solution.  Through  a  combination  of  ingenuity  and  good  fortune,  it  frequent¬ 
ly  happens  that  one  can  find  a  function  which  incorporates  the  synthetic  varia¬ 
ble  previously  derived  in  such  a  way  as  to  preserve  the  property  of  relevance 
and  yet  eliminate  the  unknown  parameters,  without  introducing  ary  new  ones.  The 
probability  distribution  of  the  resulting  quantity  will  then  bo  determinate.  A 
function  of  this  sort  is  called  a  parameter-free  variate  and  is  the  only  type 
suitable  for  the  exact  testing  cf  statistical  hypotheses  where  the  parent  distri¬ 
bution  has  unknown  parameters.  In  the  construction  of  a  parameter-free  variate, 
the  statistical  independence  of  the  component  parts  is  so  nearly  essential  that 
one  seldom  encounters  a  solvable  problem  wherein  the  components  are  not  indepen¬ 
dent.  To  derive  the  probability  distribution  of  a  function  of  two  or  more  var¬ 
iables,  their  joint  distribution  must  be  utilized  eithor  explicitly  or  implic¬ 
itly.  Now,  the  joint  distribution  of  statistically  independent  variates  can  be 
obtained  immediately  from  their  separate  distributions.  Otherwise,  the  joint 
distribution  function  would  have  to  be  known  in  its  own  right,  and  this  would  be 
a  remarkable  coincidence.  * 

Quite  understandably,  naiy  people  distrust  conclusions  based  upon  sta¬ 
tistical  evidence.  Their  suspicion  is  due  in  part  to  previous  experience  with 
fallacious  arguments  and  faulty  handling  of  data,  but  it  is  due  in  large  measure 
also  to  the  inherent  subtlety  which  characterizes  the  field  of  probability. 
Without  being  able  to  place  the  tag  of  sophistiy  upon  any  particular  step  of  a 
demonstration,  one  frequently  focls  an  intuitive  awareness  that  a  proof  is  un¬ 
sound.  Barring  act3  of  deliberate  misrepresentation,  such  a3  culling  the  data 
to  find  support  for  a  preconceived  thesis,  one  still  can  vitiate  a  statistical 
argument  by  failing  to  observe  the  conditions  upon  which  validity  depends.  But 
if  disturbing  influences  have  been  eliminated  from  the  analysis,  so  that  the 
test  variate  measures  what  it  purports  to  measure,-  if  the  correct  probability 
distribution  ha3  been  employed;  and  if  the  sample  has  been  drawn  with  due  re¬ 
gard  for  randomization — then  it  is  unreasonable  to  deny  the  conclusion  reached. 
For  under  these  conditions,  the  result  for  the  data  in  question  will  be  unique, 
and  all  competent  statisticians  will  obtain  the  same  answer. 

In  common  with  other  meteorological  elements,  rain  is  distributed  in 
a  nonrandom  fashion  both  in  space  and  in  time.  True,  the  climatological  dis¬ 
tribution  of  rainfall  at  a  given  point  Incorporates  the  existing  serial  corrc- 


lation,  so  that  the  empirical  probability  of  the  occurrence  of  ary  stated  range 
of  values  is  faithfully  represented  by  the  corresponding  climatological  rela¬ 
tive  frequency.  Nevertheless,  the  joint  probability  of  the  occurrence  of  a 
particular  set  of  values  cannot  be  derived  from  the  climatological  distribution 
of  individual  observations,  because  the  data  are  not  statistically  independent, 
end  at  the  present  tine  too  little  is  known  about  their  stochastic  relation¬ 
ships.  Therefore,  the  problem  of  obtaining  a  test  variate  having  a  known  pro¬ 
bability  distribution  is  complicated.  Vihere  individual  values  of  a  sequence 
are  not  collectively  at  random,  striking  phenomena  are  apt  to  happen.  Thus, 
apparent  periodicities  are  easily  found  in  restricted  portions  of  a  time  scries, 
and  well  defined  patterns  show  up  in  spatial  distributions.  Such  events  are 
definitely  not  in  h armory  with  statistical  independence,  and,  of  course,  a  sta¬ 
tistical  test  involving  the  tacit  assumption  of  independence  will  indicate  a 
significant  divergence  from  the  null  hypothesis.  Uith  reference  to  precipita¬ 
tion,  for  instance,  relatively  large  amounts  of  rainfall  almost  always  occur  in 
some  locality  or  other  during  nearly  every  rainstorm  of  wide  extent,  and  unless 
one  pinpoints  in  advanco  the  area  being  considered,  a  supposed  effect  can  usual¬ 
ly  be  found. 

Besides  the  difficulties  caused  by  statistical  dependence,  the  question 
of  proper  allocation  of  effects  is  a  serious  one.  Unless  the  anount  of  rainfall 
that  would  occur  naturally  can  be  taken  into  account  in  some  equitable  wy,  the 
effect  of  seeding  cannot  be  appraised. 

In  designing  experiments  to  study  the  effect  of  seeding  upon  rainfall, 
one  is  first  led  to  the  idea  of  setting  up  a  control  area  and  an  area  to  be  seed¬ 
ed.  It  would  be  advisable,  of  course,  to  determine  from  climatology  the  charac¬ 
teristics  of  these  two  areas  and  try  to  adjust  then  so  that  they  are  as  nearly 
as  possible  equally  affected  by  ary  given  rainstorm.  This  involves  the  problem 
of  how  large  the  areas  should  be  and  also  hew  close  together  they  should  be.  It 
is  often  suggested  that  one  area  be  used  as  a  prediction  for  the  second  area,  so 
that  an  estimate  night  be  obtained  as  to  what  should  occur  in  the  second  area, 
assuming  no  effect  of  the  seeding.  From  the  basic  climatology  of  the  areas  in¬ 
volved,  it  is  possible  to  compute  linear  or  nonlinear  regression  functions  con¬ 
necting  various  groups  cf  stations,  as  well  as  the  conditional  and  simultaneous 
probabilities  of  the  occurrence  of  rainfall  at  ary  two  stations  or  areas  as  func¬ 
tions  of  the  distance  between  them.  This  form  of  analysis  was  examined  in  groat 


detail  and  rejected  on  the  ground  of  one  serious  objection  which  nakes  this 
typo  of  experimentation  undesirable,  even  who n  tho  areas  to  bo  soedod  arc  docid- 
ed  each  time  by  the  toss  of  a  coin.  Although  statistically  it  is  possible  to 
compute  overall  correlations  or  probabilities  from  extended  records,  the  overall 
values  are  poor  approximations  to  current  stochastic  relationships.  Both  the 
correlations  and  the  probabilities  vary  vddely  from  year  to  year  as  functions  of 
the  i7cathcr  processes  which  arc  prevalent  at  the  tine  considered.  Consequently 
any  experinentation  utilizing  this  sort  of  technique  must  be  continued  over  a 
good  nary  years,  in  order  to  take  into  account  the  variations  in  the  parameters 
from  year  to  year.  It  nust  be  remembered  that  there  exist  such  things  as  weather 
processes  which  arc  the  dominant  feature  in  determining  tho  regression  linc3  for 
a  specific  period,  and  therefore,  the  amount  of  rainfall  at  various  stations  de¬ 
pends  upon  these  characteristics.  For  this  reason,  tho  area  being  considered 
mist  itself  act  as  a  control  area  as  well  as  an  area  of  seeding,  in  order  to  ear¬ 
ly  out  the  experiments  in  a  reasonable  length  of  tine. 

When  the  Geophysical  Research  Directorate  approached  us  with  a  specifi¬ 
cation  for  this  report,  they  intinated  that  the  work  should  be  done  so  that  the 
possibility  of  conducting  these  experiments  during  the  latter  part  of  the  sim¬ 
mer  of  1951  was  not  excluded.  For  this  reason,  the  actual  computation  and  anal¬ 
yses  were  made  on  37  years  of  data  for  the  suiner  months,  June,  July,  and  August, 
although  our  experience,  which  has  been  relatively  great,  would  indicate  that 
without  doubt  the  same  difficulties  would  occur  during  any  period  of  the  year, 
and  the  data  would  perform  analogously.  In  order  to  have  available  for  analysis 
large  areas  which  were  relatively  free  of  orographic al  and  sea  effects,  we  chose 
Iowa  and  South  Dakota  as  locations.  life  examined  very  closely  the  distribution  of 
rainfall  over  a  square  area  16$  miles  on  a  side,  whore  tho  recording  stations 
were  very  dense.  The  area  was  subdivided  into  a  square  grid  15  miles  on  a  side, 
and  the  rainfall  values  at  the  intersection  points  of  the  grid  were  estimated  by 
interpolation.  This  was  done  for  hourly,  weekly,  an£  monthly  observations.  The 
emphasis  was  placed  upon  hourly  data,  since  tho  behavior  of  individual  storms 
could  be  studied  most  effectively. 

Attempts  woro  made  to  characterize  the  precipitation  patterns  by  util¬ 
izing  orthogonal  polynomials — a  technique  which  had  proved  fruitful  in  previous 
work  with  pressure  maps.  Unfortunately,  the  precipitation  patterns  wore  so  com¬ 
plicated  that  a  workablo  number  of  polynomials  failod  even  to  represent  more  than 
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70  percent  of  the  variability,  even  when  the  data  vrcre  smoothed  in  a  reasonable 
fashion  .  Furthcmorc,  when  largo  /mounts  of  rainfall  were  added  to  various  sec¬ 
tions  of  the  grid,  in  order  to  represent  a  lypothetical  seeding  effect,  it  was 
not  possible  to  observe  a  characteristic  change  in  the  polynomials.  Consequent¬ 
ly,  the  polynomials  could  not  bo  usod  to  indicato  changes  in  the  rainfall  dis¬ 
tribution. 


It  would  readily  occur  to  anyone  that  tad r king  with  smaller  areas  would 
mitigate  the  difficulty  of  characterization,  and  that  perhaps  simple  tests  of 
significance  could  then  be  devised.  Granting  this,  it  is  still  difficult  to 
determine  what  size  the  area  should  be;  for  there  must  be  room  enough  to  3eed  one 
portion  and  reserve  the  rest  for  control,  and  yet  the  area  must  be  compact  enough 
to  permit  a  feasible  analytical  representation  of  the  observed  rainfall,  under 
natural  conditions.  However,  it  is  along  this  direction  that  we  finally  decided 
upon  a  test  which  appears  to  be  satisfactory.  Standard  experimental  designs  were 
critically  examined,  and  even  some  nonlinear  hypotheses  were  explored.  Particu¬ 
lar  consideration  was  given  to  these  designs  which  lend  themselves  to  the  use  of 
mean  square  successive  differences  and  other  schemes  for  eliminating  pronounced 
local  trends  in  the  rainfall  pattern.  All  of  these  experimental  designs  were 
found  to  be  unsound  in  the  light  of  the  observed  behavior  of  rainfall,  if  the 
area  were  at  all  large. 

After  the  examination  of  a  great  doal  of  data  in  the  area  chosen,  the 
plan  presented  in  the  Appendix  appeared  to  be  the  only  workable  one  which  com¬ 
bined  a  tenable  hypothesis  with  a  sufficient  reduction  of  the  meteorological 
variability,  so  that  if  seeding  has  arc-  effect  it  should  be  observable  with  a 
reasonable  number  of  experiments.  Particular  details  of  the  experiment  are  open 
to  change  and  can  be  modified  in  v/eys  that  seem  practical.  It  is,  of  course,  pos¬ 
sible  that  the  actual  seeding  cannot  be  performed  adequately  for  this  type  of  con¬ 
figuration  and  that  adjustments  will  have  to  be  made  in’ the  light  of  practicality 
which  will  reduce  the  efficiency  of  the  teat. 


In  the  method  ultimately  developed,  smoothing  was  not  employed. 
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Thc  final  plan,  which  seems  most  adaptable  to  the  experiments,  utilizes 
a  square  grid  containing  nine  observation  stations,  located  about  IS  miles  apart. 
The  station  at  the  center  of  the  grid  and  the  one  in  the  middle  of  the  edge  which 
is  directly  downwind  arc  the  two  stations  to  be  seeded.  The  actual  rainfall  to 
bo  expected  at  those  two  stations,  exclusive  of  seeding  effects,  can  be  estimat¬ 
ed  very  well  from  the  remaining  seven  stations.  Therefore,  the  deviations  of 
tho  actual  from  estimated  values  at  these  two  stations  provide  the  basis  for  the 
test  variate. 

In  our  view,  this  design  eliminates  the  cause  of  statistical  dependence. 
'.Te  identify  the  spatial  correlation  from  point  to  point  with  an  underlying  con¬ 
tinuous  geometrical  distribution  of  rainfall  existing  at  the  time  of  observation, 
and  we  attribute  the  serial  correlation  to  the  pcr3istenco  of  this  pattern. 
Therefore  wc  have  hypothesized  that  the  removal  of  the  geometrical  idealization 
would  leave  residuals  which  are  individually  and  collectively  at  random.  Lack¬ 
ing  rainfall  observations  located  exactly  in  a  square  array,  we  chocked  tho  hy¬ 
pe  ti  isized  distribution  of  the  test  variate  as  best  we  could  with  data  read  from 
contour  maps  of  actual  rainfall.  !r7e  wish  to  call  explicit  attention  to  this  fact 
and  recommend  that  a  few  trials  be  made  with  direct  observations  without  seeding. 
Using  $0  sets  of  nine  values,  we  found  that  the  Kolmogorov-Smimov  te3t  for  good¬ 
ness  of  fit*  supported  tho  hypothesis  at  better  than  the  20  percent  level  of  sig¬ 
nificance. 


Uaty  of  these  sets  were  taken  from  contiguous  (though  nonoverlapping) 
areas  and  from  consecutive  hours,  but  we  advise  a  different  sampling  procedure 
in  field  practice.  If  a  great  many  observations  vroro  taken  over  the  experimen¬ 
tal  area,  the  shape  of  the  geometrical  surface  could  then  be  determined  within 
any  preassigned  tolerance.  In  that  event,  it  should  not  matter  how  dense  the 
samples  are  in  spaco  or  time.  As  a  practical  expedient,  however,  wo  have  approx¬ 
imated  the  true  surface  by  a  properly  oriented  plane.  This  approximation  should 
be  quite  satisfactory  for  the  limited  area  considered,  but  it  introduces  the  pos¬ 
sibility  of  a  small  systematic  error,  if  the  samples  are  too  close  together. 
Therefore,  randomization  should  be  employed  in  actual  operations. 


*Frank  J.  Massey,  Jr. »  "The  Kolmogorov-Smimov  Test  for  Goodness  of  Fit", 
Journal  of  the  American  Statistical  Association,  Vol.  U6,  No.  253,  Uarch,  1951, 
p.  68-78. 


With  regard  to  tho  degree  of  approximation  attainable  with  a  plane  sur¬ 
face,  we  have  found  that  the  multiple  correlation  is  very  high.  In  about  35  per¬ 
cent  of  the  cases,  the  multiple  correlation  exceeded  . 95;  in  about  55  percent,  it 
oxcccdcd  .90;  in  about  65  percent,  it  exceeded  .85;  and  in  about  75  percent,  it 
exceeded  ,80.  Accordingly,  it  would  be  feasible  to  restrict  tho  formal  analysis 
to  those  samples  for  which  the  plane  surface  is  a  good  fit.  This  can  bo  done  by 
having  a  fow  samples  to  spare*  and  it  does  not  violate  any  principles  of  rigor, 
because  the  selection  will  be  made  entirely  without  reference  to  the  amounts  of 
rainfall  at  the  test  points. 

If  seeding  increases  rainfall  to  any  appreciable  extent,  the  test  vari¬ 
ate  should  respond  unmistakably.  The  kind  of  rosponso  available  in  a  typical 
sample  is  exhibited  in  Table  I.  Here  the  value  of  the  test  variate  was  computed 
in  five  cases  from  rainfall  data  obtained  as  previously  stated. 

Table  I 


Values  of  Test  Variate  Under  Three  Conditions 


Case 

Natural 

Conditions 

Addition  of  ,Q5  in. 
at  Test  Points 

Addition  of  .10  in. 
at  Test  Points 

1 

CM 

• 

1 

3-UO 

6.17 

2 

1-75 

2.90 

3.92 

3 

-.9U 

U.91 

9.69 

h 

.U5 

3.33 

5.U2 

5 

-.76 

7.83 

13.60 

Then  the  original  rainfall  at  the  two  hypothetically  seeded  points  was  arbitrar¬ 
ily  increased  by  .05  in.,  and  the  corresponding  value  of  the  test  variate  was  re¬ 
computed;  finally,  the  original  rainfall  at  the  two  points  was  arbitrarily  in¬ 
creased  by  .10  in.,  and  tho  corresponding  value  of  the  te3t  variate  was  conputod 
again.  The  test  variate  is  strongly  affected  in  every  case,  although  not  to  a 
uniform  degree,  inasmuch  as  the  internal  variability  of  the  actual  data  is  also 
involved.  The  probability  levels  associated  with  each  value  of  the  test  variate 
arc  shown  in  Tabic  II,  Even  the  individual  values  obtained  from  the  modified  rain- 
fall  amounts  are  significant;  the  collective  significance  levels  ,  however,  would 

*The  collective  significance  level  .is  not  the  product  of  the  separate  probabili¬ 
ties,  For  exposition,  see  Appendix, 
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Table  II 

Associated  Probability  Levels 


Case 

1 

.65 

.oil* 

.0018 

2 

.08 

.022 

.0087 

3 

O 

00 

• 

.001* 

.0003 

1* 

.31 

.oil* 

.0028 

5 

.75 

.001 

.0001 

be  utterly  conclusive.  For  natural  conditions,  the  collective  significance  level 
is  .52,  which  comfortably  sustains  the  null  hypothesis;  for  the  addition  of  .05  in. 
the  collective  significance  level  is  about  2.1*  x  lCT?,  which  indicates  that  the 
null  hypothesis  is  most  definitely  untenable;  and  for  the  addition  of  .10  in., 
the  collective  significance  level  is  about  8.5  x  10"^-,  Trtiich  is  so  small  that  it 
•would  be  preposterous  to  entertain  the  null  hypothesis, 

Che  night  ask  whether  .05  in.  and  .10  in.  represent  uncommonly  large 
hourly  rainfall  amounts  in  themselves.  In  the  area  considered,  the  empirical 
distributions  of  hourly  rainfall  during  the  month  of  July  were  obtained  for  five 
regular  weather  stations.  Table  III  presents  tho  relative  frequencies  with  which 
the  hourly  rainfall,  when  it  docs  rain,  equals  or  exceeds  ,05  in.  and  .10  in.  re- 

Table  III 

Relative  frequencies  with  which  July  Rainstorms  at  Five  Stations  in  Iowa 
Produce  Hourly  Amounts  of  Stated  Magnitudes 


Station 

At  Least 
.05  in. 

At  Least 
.10  in. 

At  Least 
.15  in. 

At  Least 
.20  in. 

Burlington 

sT " 

"’'33 

.26 

.21 

Coon  Rapids 

.u* 

.26 

.18 

.13 

Des  Moines 

.U3 

.29 

.20 

.11* 

Keokuk 

.U9 

.33 

.20 

.16 

Washington 

.56 

.37 

.26 

.19 

spectively,  as  Troll  as  the  sane  information  for  .15  in.  and  .20  in.  from  tho 
table,  it  is  evident  that  an  hour's  accumulation  of  rain  frequently  amounts  to  as 
much  as  .05  in.  or  .10  in* 
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Becauso  tho  proposed  experimental  areas  are  sufficiently  small,  they 
can  be  located  in  many  places,  and  ary  number  of  them  can  be  pooled  to  obtain  a 
combined  measure  of  significance,  as  described  in  the  Appendix,  Although  this 
set-up  is  probably  not  the  only  one  >7hich  night  pork,  nc  hold  that  any  valid  sys¬ 
tem  would  have  to  preserve  tho  essential  attributes  of  this  design. 


The  contemplated  statistical  test  calls  for  a  criterion  variate  which 
is  approximately  normally  distributed.  Hourly  rainfall  in  a  cross  between  a,  dis¬ 
crete  and  continuous  variate^  having  a  fairly  large  probability  of  the  occurrence 
of  exactly  zero,  an  appreciable  probability  of  a  trace  of  rainfall,  and  in  the 
measurable  range,  a  skewed  probability  density  function  of  rainfall  as  a  contin¬ 
uous  variate.  To  obtain  a  statistically  manageable  variate,  we  shall  exclude 
from  consideration  values  of  hourly  rainfall  less  than  . 0 $  in.  and  to  reduce  the 
skewness  we  shall  work  with  the  natural  logarithm  of  the  observed  rainfall. 


As  a  compromise  between  the  requirements  of  keeping  the  test  area  small 
enough  to  be  statistically  homogeneous  and  large  enough  to  support  a  subdivision 
into  seeded  and  nonseeded  portions,  we  have  chosen  to  vrork  with  experimental  plots 
containing  nine  observation  points.  Ideally  these  plots  should  be  square  in  shape 
and  subdivided  into  quadrants,  with  the  observations  lying  at  the  corners  of  the 
grid  squares,  as  shown  in  Figure  1.  The  actual  orientation  of  the  plot  must  be 


In  Figure  1,  the  symbols  (i,j  *  1,  2,  3)  represent  tho  nine  observa¬ 
tions  of  rainfall  and  indicate;  at  the  same  time;  the  locations  of  tho  observa¬ 
tion  points.  Choose  a  rectangular  coordinate  system  with  origin  at  the  center  of 
tha  square  (R^),  tho  X-axia  coinciding  with  the  lino  R?1R22R23’  the  f-axia 
coinciding  with  the  line  ^2^22^32*  ^a^nS  the  side  of  a  grid  square  as  the  unit 
of  length,  wb  set  up  the  X-Y  coordinate  pattern  shown  in  Figure  2. 


Jrilil _ iSOl _ 0*1) 
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Figure  2 


Within  the  test  plot,  the  natural  logarithm  of  the  rainfall  at  the  point 
i»  assumed  to  satisfy  a  linear  regression  equation  given  fcy 

+i}  •  K  *  ca  * 

Here  wa  make  the  standard  assumptions,  namely  that 

rij  5  pij  +  Hi 

where  the  *»s  are  Independent  normally  distributed  variates  with  zero  mean  and 
2 

cannon  variance  a  . 

i 

Throe  points  should  be  borne  in  mind  with  regard  to  the  regression  func¬ 
tion  p.  First,  end.  most  important,  lm  the  feet  that  we  do  not  have  to  assume  that 
the  parameters  p^  Pj,  pg  are  invariant  from  place  to  place  on  the  nap,  or  from 
storm  to  storm  at  a  fixed  location*  Whatever  the  parameters  might  be  for  the 
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storm  and  test  plot  in  question,  to  shall  estimate  them  from  the  observational 
data  taken  from  that  one  time  and  place,  vrithout  dependence  upon  climatological 
records.  Second,  the  statistical  test  of  significance  will  be  parameter-free,* 
that  is,  it  will  bo  a  determinate  function  of  observable  quantities.  Third, 
even  if  the  distribution  of  the  logarithm  of  rainfall  is  not  very  nearly  normal, 
it  is  still  quite  reasonable  to  assume  a  normal  distribution  for  the  residual 
variate  z,  because  the  removal  of  a  definite  geometric  pattern  p  from  r  could 
well  leave  normally  distributed  random  deviations. 

The  hypothetical  regression  function  p  will  be  approximated  by  a  sim- 
ilar  function  r  where 


r*  "  bo  +  *1*  +  V 

and  the  constants  bQ,  bj,  b2  will  be  determined  from  the  data.  Seeding  will  be 
confined  to  two  points,  (0,0)  and  (0,1).  The  remaining  seven  points  will  be  used 
to  fit  the  constants.  Since  these  points  will  not  be  seeded,  the  function  r 
will  provide  an  estimate  of  the  amount  of  rainfall  which  would  occur  naturally. 
Briefly,  the  object  of  the  dosign  is  this :  The  experimental  control  will  be  furn¬ 
ished  by  the  regression  function,  and  the  measure  of  effectiveness  will  be  based 
upon  the  deviations  of  actual  rainfall  from  the  regression  function  at  the  two 
seeded  points.  The  experiment  will  be  replicated  as  often  as  necessary  to  build 
up  a  conclusive  body  of  evidence.  The  statistical  theoiy  and  computational  pro¬ 
cedure  are  presented  in  the  following  paragraphs. 


From  established  principles  of  mathematical  expectation,  it  can  be 
proved  that  an  unbiased  estimate  of  the  residual  variance  o^  is  given  by  s^t 
where 


s 


2 


2 

Being  a  function  of  random  variables,  the  statistic  3  ,is  itself  a  random  variable. 
Its  probability  distribution  is  simply  related  to  that  of  chi-square,  inasmuch  as 
tho  variate  Ub^/c^  is  distributed  precisely  as  chi-square  with  four  degrees  of 
freedom. 


If,  under  the  assumptions  stated,  a  statistically  derived  regression 

M 

function  r  is  applied  to  fresh  data,  the  attendant  deviation  -  r  jj,  regard¬ 
ed  as  a  random  variable,  will  be  normally  distributed  with  zoro  mean,  but  its 
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variance  will  be  a  function  of  the  sampling  variances  and  covariances  of  the  es¬ 
timated  regression  constants  b^  b^,  bg.  The  variance  of  -  r*^j  is 


E(riJ  "  r 


id 


r  -  E(r. 


id 


-  'id>‘ 


+  E(r 


id 


-  Oij)' 


"*  *  «r*«  -  '1/ 


5y  direct  substitution 


e(p*lj  -  l’i})2  *  «<».  -  V  ♦  (h  -  h*±  *  <b2  -  h*}? 

To  evaluate  this  Quadratic,  matrix  notation  is  convenient*  Introduce  the  matrix 


N 

2x 

Sy 

A  - 

2x 

Sc2 

zxy 

2* y 

5r2 

and  denote  its  inverse  by  A”\  The  first  element  N  of  the  matrix  A  represents 
the  sample  size,  -which  in  this  case  is  7«  Define  the  row  vector 

v  -  (1  x±  yj) 

and  let  v*  stand  for  its  transpose.  Then  it  can  be  demonstrated  that 

«(b#  -  PQ)  *  -  Pj^  ♦  0>2  -  P2^^2  "  0,2  ▼A"**1 

Hence  the  variance  of  r^-r*^  is 

E(pid  **  r*id)2  "  02(1  +  vArly,) 

The  latter  quantity  in  parentheses  is  a  numerical  constant  which  can  be  computed 
without  reference  to  the  rainfall  amounts. 

Consider  now  the  deviations 

“o  =  r22  -  r* 22 
\  ■  r12  "  r<12 


2  '/ 

The  variance  of  u  is  k  a  and  that  of  u.  is  1^,  cr  ‘  where 

0  0  1  J. 

k  -  (1  +  v  A_1v*  ) 
o'  o  o 

*  (1  +  v^A'^v'^) 

in  which  vq  -  (1  0  0),  ■  (1  0  1).  Vfe  have  observed  that  uq  and 

are  normal  variates,  each  having  a  mean  of  zero.  If  wn  now  construct  the  linear 

combination 

.  _  uo//% + 


this  too  will  be  normally  distributed  with  zero  mean  and  its  variance  will  be 


cr2/  2. 


The  two  statistics  and  u  are  independently  distributed.  Accordingly, 


if  we  set 


yy° 72  3  v? 

y  s2/o2  y ? 


we  arrive  at  a  parameter-free  statistic  having  the  well-known  t-distribution  with 
four  degrees  of  freedom.  In  other  words,  the  probability  density  function  of  t 


,  -  -5/2  2-5/2 

f(t)  -  |(1  ♦  tzA)  -  12(1*  +  tz) 


If  seeding  increases  rainfall  then  u,  and  consequently  t,  should  tend  to  exceed 
zero.  Therefore,  the  appropriate  measure  of  significance  is  the  probability 
(under  the  null  hypothesis)  that  a  value  of  t  at  least  as  great  algebraically  as 
the  one  obtained  in  the  seeding  experiment  would  occur  by  chance.  Specifically, 
if  in  the  experiment  we  find  that  t  ■  x,  then  trie  measure  of  significance  p(*r)  is 
given  by 


p(v)  •  700 f(t)dt  ■  ^  [1  -  (.fo.  -■  y \t>g-  ] 

(1*  ♦  %  ) 


While  the  foregoing  equation  gives  the  exact  value  of  p,  numerical  work  can  be 


& 

reduced  by  referring  to  suitable  tables.  The  curve  of  Figure  3  represents  the 
density  function  f(t),  and  the  shaded  portion  under  the  curve  represents  the  prob¬ 
ability  p(*t). 


Figure  3 

The  experiment  can  be  extended  indefinitely  by  replication — that  is, 
applying  the  same  scheme  over  and  over.  Several  nonoverlapping  plots  can  be  used 
in  the  same  storm,  and  repeated  samples  can  be  taken  from  the  same  plot  in  differ¬ 
ent  storms.  The  total  evidence  can  then  be  pooled  by  the  standard  procedure  of 
combining  probabilities  from  independent  tests.  Although  the  technique  is  simple, 
and  to  a  statistician  intuitively  clear,  the  non-statistician  is  apt  to  experience 
difficulty  in  appreciating  its  validity.  Therefore  we  shall  sketch  a  derivation  of 
the  method  from  first  principles. 

Suppose  that  tvro  independent  tests  resulted  in  the  significance  measures 
(probabilities)  p^  and  p2  respectively.  Loosely  speaking,  the  probability  of  ob¬ 
taining  such  an  outcome  by  chance  is  the  product  P^P2*  tore  precisely,  however, 
this  product  is  really  a  probability  density,  because  and  p2  are  both  continu¬ 
ous  variates.  As  a  matter  of  fact,  it  can  be  shown  that  and  p2  are  rectangular¬ 
ly  distributed,  and  since  they  obviously  range  from  zero  to  1,  the  constant  density 
value  of  each  is  unity.  Construct  a  rectantular  coordinate  system,  as  in  Figure  k, 
with  horizontal  axis  p^  and  vertical  axis  p2#  Under  the  null  hypothesis  the  points 
(Pl,P2),  representing  the  results  of  any  pair  of  tests,  pre  uniformly  distributed 
over  the  unit  square  in  the  first  quadrant.  The  rectangular  hyperbola  p^p2  ■  X 
defines  the  probability  locus  of  all  pairs  of  test  results  having  equal  likelihood. 
The  combined  significance  measure  of  any  pair  of  tests  for  which  p^p2  -  X  is  equal 
to  the  area  of  that  portion  of  the  unit  square  which  lies  below  the  hyperbola,  for 

*1L  G,  Kendall,  The  Advanced  Theory  of  Statistics,  Vol,  1;  London,  Charles  Griffin 
and  Co»  Appendix,  Tabic  3T  p.LtUO,  The  numbers  in  this  table,  for  four  degrees  of 
freedom,  give  the  appropriate  significance  probabilities  when  subtracted  from 
unity* 


Figure  U 


this  area  represents  the  total  probability  of  obtaining  some  pair  of  values  for 
which  the  likelihood  is  equal  to  or  less  than  X.  Denoting  the  combined  signifi¬ 
cance  by  P,  ire  have  ^ _ 

P  m  1  -  £  J\/p^  dPidP2  "  M1  “  lnX  ) 

Similarly,  for  three  independent  tests,  the  sample  space  is  a  unit  cube,  within 
which  all  points  have  unit  density.  The  likelihood  contour  is  the  hyperbolic  sur¬ 
face  p^PgP^  ■  X,  and  the  combined  significance  measure  is 

p  • 1  *  t  -  x  [i  -  mx. 

In  general,  for  n  tests  it  is  easy  to  show  that  the  combined  significance  measure 

P  -  X(1  -  in X  +  &$£  _  ♦  ...  +  (-if-1 

* 

where  X  equals  the  product  of  the  p’s.  Since  very  good  tables  of  logarithms  are 
available,  P  can  be  computed  directly  from  the  foregoing  equation.  However,  most 
of  the  arithmetic  can  be  avoided  by  using  tables  of  the  chi-square  integral.  For 
2n  degrees  of  freedom,  the  density  function  of  chi-square  is 


g(x2) 


1 

2°r(n) 


(tr1  e-*2/* 


and  it  turns  out  that  the  corresponding  probability  integral  is 


jf  g(x2)<fc2  -  •"*  (1  *  (x2/2)  +  *  ^ . gftH 


,2  In  \1*-1 


-y^/2 

Now,  if  we  make  the  Substitution  X  a  e  A  '  we  obtain 


P  .  e-^tl  ♦  (x2/P)  *  ♦  gjpt  ♦...*$$}  1  *  n  8(X2)  dx! 


Therefore,  the  combined  significance  measure  P  can  be  evaluated  from  the  integral 
of  chi-square  with  2n  degrees  of  freedom  if  we  define 

X2  -  -21n  X 

A  statistician  arrives  at  this  conclusion  much  more  expeditiously  by  associating 
with  each  probability  p,  a  chi-square  equal  to  -21n  p.  and  having  two  degrees  of 
freedom,  by  the  reproductive  property  of  chi-square,  the  sum  of  these  is  another 
chi-square  but  with  2n  degrees  of  freedom.  The  step  of  defining  an  individual 
chi-square  as  -21n  is  justified  by  the  fact  that  for  two  degrees  of  freedom, 
chi-square  1 a  precisely  equal  to  minus  twice  the  logarithm  of  its  own  probability 
integral* 

Proceeding  now  to  specific  computational  details,  we  present  a  comprehen¬ 
sive  numerical  routine  for  tho  proposed  grid  system.  The  corresponding  values  of 
Xj,  y^,  and  r^  arc  exhibited  in  Table  1,  and  the  entries  of  A  are 

M  ■  7#  2x  •  0,  ?y  -  -1,  Sc2  ■  6,  2xy  -  0,  by2  ■  5 


Bence 


7  0-1 

0  6  0 

-10  5 


and  the  inverse  matrix  is 


_  1 
A  “ss; 


30  o  6 

0  3U  0 

6  o  1*2 


/f-?- 


Tablo  1 

Corresponding  Values  of  Variates 


iT~ 

xi 

y'j  4y 

li 

-1 

1  rn 

13 

1 

1  r1? 

21 

-1 

0  r21 

23 

mL 

0  v 

r2.3 

31 

-1 

r 

r3l 

32 

0 

-1  r32 

33 

1 

-i  r33 

The  matrices  A  and  A“^  are,  of  course,  independent  of  the  rainfall  data, 
tion,  wo  need  four  functions  of  the  observations — namely: 

co  ■  ■  <rll  *  ra  *  r31)  *  (r13  *  r23  +  r33>  *  r32 

°1  ’  airiJ  •  -<ru  *  r21  *  r31>  *  (r13  *  r23  *  r33> 


c2  *  »Jril  '  (rU  *  r13>  -  <r31  *  r32  *  r33> 


2222222! 

*13 .  “  rn  .  +  r21  +  r31  *  r13  +  r23  +  r?3  +  r32 

Formally,  the  estimated  regression  constants  arc  given  ly 


In  addi 


”bo~1 

V 

*1 

.•1 

«  A 

C- 

.i. 

b2 

°2 

M  Wi 

•  mm 

-co  +  c: 

— 3IT* 


In  this  case  the  solution  is 


A  -11- 


The  remaining  quantities  are 

.  vA  *  vA  _  s/*5 « 

u  2  2 - 

and 

t.iZL 

JJ 

_  2 

If  no  interest  attaches  to  u  and  a  in  themselves,  the  calculation  of  t  can  be 
shortened  by  combining  terms  in  a  different  wqy  obtaining 

w  u  ♦  w,u, 

t  -  -  p  0  ■  3-7  - 

y  Z( ru  -  r\~)2 

where  _  _ 

wo  "  J2\  "  v/68/39  *  i*32 

W1  “  M  '  y^U  -  1.19 

M  p  O 

and,  of  course,  Zfr^  “  r  “  (hQc0  +  '’I®!  +  ^2C2^*  To  six  decimals  the 

weights  irQ,  w^  are  1. 3201*51*  1.190238  respectively.  Thus  the  figures  given  above 
are  actually  good  to  three  decimals,  which  should  suffice  for  practical  purposes. 

Operational  considerations  might  dictate  tho  use  of  a  fixed  network  of 
points  for  all  experiments.  In  that  case  the  test  plots  could  be  chosen  as  small 
compact  areas  containing  nine  points,  with  tho  two  test  points  oriented  as  before. 
An  important  condition  to  be  satisfied  is  that  none  of  the  nonscedod  points  must 
be  downwind  from  the  test  points.  Because  the  network  is  no  longer  assumed  to  be 
strictly  regular,  the  double  subscript  notation  is  inappropriate.  Instead,  wb 
shall  number  the  points  0,  1,  2, 8,  the  first  two  being  the  test  points.  This 
scheme  is  represented  in  Table  2. 
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Table  2 

fixed  Network  Scheme 


Serial 

Number 

X 

7 

r 

0 

xo 

y0 

r7 

1 

*1 

yi 

rl 

2 

X2 

y2 

*7 

3 

X3 

y3 

r3 

h 

xh 

yU 

rl; 

5 

*5 

y5 

r5 

6 

*6 

y6 

r6 

7 

*7 

y7 

r7 

8 

x8 

y8 

r8 

Seeded 


Nonaeeded 


Although  the  regression  constants  b^,  b^,  b^  and  the  inverse  matrix  aw  formal- 
involved  in  the  determination  of  s2,  u,  and  t,  they  need  not  be  derived  explicit- 


•  An  efficient  computational  procedure  is  as  follows: 

Compute  the  sums  of  first  powers,  squares,  and  cross  products — 


8*8  8 


8 


8 


8  p  8  9  ®  2  ® 

f*i*  2^1'  f'v  2^i’  *i  *  ^*iyi'  ^iri»  ^iri 


Then  set  up  the  initial  matrix  as  shown  bolcw,  it  being  understood  that  all  sum- 
nations  go  from  2  to  8, 


7 

Zx 

Zy 

Zr 

1 

.  1 

2x 

Zx2 

ay 

Zxr 

*o 

X1 

Initial 

2Sr 

Zxy 

3 r2 

Zfcrr 

yo 

yl 

Matrix 

Zr 

Zxr 

Zjrr 

— 

_ 

— 

In  the  loner  right  hand  corner  of  the  initial  matrix,  the  dashes  represent  blanks, 
which  are  to  be  ignored  in  the  computations.  Now  derive  the  Crout  auxiliary  mat¬ 
rix  according  to  the  directions  given  in  Marchant  Methods,  MU  182.  Representing 
the  general  element  of  this  auxiliary  matrix  by  tho  symbol  we  have 


Auxiliary 

Matrix 


°11 

e12 

e13 

elli 

ei5 

el6 

C21 

°22 

e23 

C2h 

e25 

°26 

C31 

e32 

°33 

e3U 

c35 

e36 

V 

%2 

Cii3 

— 

— 

— 

e5l 

c52 

e53 

— 

— 

— 

°61 

C62 

e63 

— 

— 

— 

Then  2  8  *  2  8  2 

Us  -  -  r  A)  -  jr±  -  (c^e^  +  e^  ♦  e^e^) 

uo  "  ro  -  r*o  ’  ro  -  <c5?lell*  +  e52e2U  +  e53e3U} 

Ui  "  rx  -  r\  -  rx  -  (o6le^  +  e^e^  ♦  e^e^) 

ko  “  1  +  c5iei5  +  e52e25  +  e53B35 

*1  "  1  +  e6lel6  +  c62e26  +  e63e36 

The  rest  of  the  calculations  arc  performed  as  previously.  The  mathematical  justi¬ 
fication  of  tliis  computational  short-cut  requires  too  much  matrix  theory  to  be 
gone  into  here,  but  a  numerical  verification  is  furnished  by  the  following  cxarfJe; 

A  square  plot,  30  miles  on  a  side,  was  chosen  in  central  Iowa.  As  pre¬ 
viously  discussed,  this  was  subdivided  into  four  grid  squares  15  miles  on  a  side. 
From  a  contour  map  of  hourly  precipitation  as  of  U  a.m.  on  18  June  1950  the  rain¬ 
fall  amounts  were  read  at  the  intersection  points  of  the  grid.  The  actual  values 
are  shown  in  Figure  5a,  and  the  corresponding  natural  logarithms,  each  inemasod 
by  10,  are  indicated  in  Figure  5b.  It  is  legitimate  to  add  10,  or  ary  other  ar¬ 
bitrary  constant,  to  the  logarithms,  because  the  regression  function  adjusts  to 
ary  linear  change  of  variable;  one  should  be  careful,  however,  not  to  forget  to 
add  10  to  the  logarithm  when  the  rainfall  amount  exceeds  unity.  17b  shall  perform 
the  computations  first  by  tho  method  given  for  the  square  grid  and  then  by  the  gen¬ 
eral  method  using  the  Crout  auxiliary  matrix.  In  ordor  to  insure  numerical  agree¬ 
ment  between  the  two  systems,  we  shall  carry  more  figures  than  would  otherwise  be 
needed. 
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8.33 9  0.285  8.285 


Since  all  nine  of  the  observations  represent  natural  precipitation,  we  should  ex¬ 
pect  the  value  of  t  to  support  the  null  hypothesis  that  ordinary  circumstances  are 
in  force. 


The  functions  c^  Cp  etc.  are 


cQ  -  60.11*9 
cx  -  -.1*29 
c2  -  -9.721* 

Zr2  -517.11*621*1 


The  regression  constants  are 

5c  +  c 


b  - 


5J— ^  -  22^021  -  8.5591*1*1 


bl  -  ^  -.071500 

b2 "  "  -• 232912 


The  sun  of  squares  of  residuals  is 


fi  -15- 


Z(r-r*)2  -  Zr2  -  (bco  +  +  b2c2) 

-  517.1U62U1  -  517.137326 


and 


-  .008915 
ff2  -  2(r-r*)2  . 

4 


.002229 


The  deviations  uq  and  u^  and  the  linear  combination  u  are 


u  -  r99  -  b  -  8.530  -  8.559Uil  =■  -.02914*1 
o  zz  o 


^  -  r12  -  b0  —  bj  -  8.285  -  8.326529  -  -01il529 


-  ..*J*?*  «  _  -.asaig  _ 


031220 


In  connection  with  u,  we  note  for  future  reference  that 

kQ  -  39/3U  -  1.11*7059,  -  214/17  -  1.1411765 

Now  we  compute  the  value  of  t  by  the  two  formulas  given  above,  with  and 

_  0 

without  explicit  use  of  u  and  s  .  I$r  the  first  formula, 

.  _  u/2  _  -.031220  \fi~_  -.014*152  _ 
t  -  — 7CC7212 —  ".'0*732  -935Z 

I$r  the  second  formula,  using  w  ,  w^  to  six  decimals  (wQ  ■  1.3201*51*  w^  ■  1.190238) 
we  have 


yz(r  -  r*)2  -  y. 008915  -  .0914*19 


.  _  Vo  +  Vl  _  -.088305  . 

t  ■  -7 . -  -9352 

y/z(r  -  rj 


Again,  by  the  second  formula  but  using  wQ*  w^  to  two  decimals 


.  _  ;o88?8 

*  m7ms 


.9350 


Thus  the  two-decimal  weights  are  sufficiently  accurate. 


The  general  method  is  designed  to  take  care  of  irregularly  spaced  points 
but,  of  course,  it  will  give  the  same  answer  as  the  first  method  where  the  latter 
is  applicable.  The  arrangement  of  data  for  the  general  method  is  shown  in  Table  3* 


Table  3 

Basic  Data  for  General  Method 


The  corresponding  Crout  auxiliary  matrix  is 


